Outcomes and Duration of Tidal Evolution in a Star-Planet-Moon System 



Takashi Sasaki 

Department of Physics 

University of Idaho 
Moscow, ID 83844-0903 
ResearcherID:B-44 70-201 1 
t s as aki@ vandals . uidaho . edu 



(N 

o 

(N 
G 

<N 

w 

Oh 

6 



> 

m 
m 
o 

\6 
o 

(N 

> 



Jason W. Barnes 
Department of Physics 

University of Idaho 
Moscow, ID 83844-0903 
ResearcherlD: B-1284-2009 
jwbarnes@uidaho . edu 

David P. O'Brien 

Planetary Science Institute 
1700 East Fort Lowell, Suite 106 
Tucson, AZ 85719-2395 
obrienOpsi . edu 



ABSTRACT 

We formulated tidal decay lifetimes for hypothetical moons orbiting extrasolar planets with both lunar and 
stellar tides. Previous work neglected the effect of lunar tides on planet rotation, and are therefore applicable only 
to systems in which the moon's mass is much less than that of the planet. This work, in contrast, can be applied 
to the relatively large moons that might be detected around newly-discovered Neptune-mass and super-Earth 
planets. We conclude that moons are more stable when the planet/moon systems are further from the parent 
star, the planets are heavier, or the parent stars are lighter. Inclusion of lunar tides allows for significantly longer 
lifetimes for a massive moon relative to prior formulations. We expect that the semi-major axis of the planet 
hosting the first detected exomoon around a G-type star is 0.4-0.6 AU and is 0.2-0.4 AU for an M-type star. 

Subject headings: celestial mechanics-planetary systems-planet and satellites:general 



INTRODUCTION 



The first discovery of an extrasolar planet in orbit 



around a main-sequence star was made by Mayor & Queloz 
( |1995 ). Since then, more than 700 extrasolar planet^ 
have been discovered. Although extrasolar moons have 
not yet been detected, they almost certainly exist. Most 
of the planets in our solar system have satellites. Even 
Pluto, though no longer officially a planet, has three moons 
(Weaver et al. 112006 1. It is likely that the mechanisms for 



Robutel|1 993). Stable obliquity in its star's habitable zone 
may be necessary for a planet to support life. An Earth- 
size planet with no moon, or a relatively small one, may be 
subject to large fluctuations in obliquity. In such a case, 
favorable conditions may not last long enough for life to 
become established. In the same way, orbital longevity 
is required for any life form to have time to become es- 
tablished. Hence, the prospects for habitable planets may 



hinge on moons (Ward & Brownlee 20001; but see also 
( |Lissauer et"aLfl2012| ). 



moon formation in our solar system (impact, capture, and 
coaccretion) prevail beyond it( Mamajek et aL]|2011 1. 



In 2005, Rivera et al. (2005) discovered Gliese 876 d, 



The Earth's obliquity, or axial tilt, is stabilized by the 
Moon (Laskar et al. 1993). Mars, on the other hand, has 



relatively small satellites, and its obliquity changes chaot- 



thc first super-Earth around a main sequence star. To date 
more than thirty super-Earths have been discoverecj^] The 
discovery of Kepler 22-b in the habitable zone gives rise to 
the possibility of life beyond our Solar System (iBorucki 



ically, fluctuating on a 100,000-year timescale (Laskar & 



et al.||2012[ ). It is important to know the lifetime of moons 



1 http:/ /exoplanet.eu/ 



2 http:/ /exoplanet.eu/ 
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orbiting super-Earths in the habitable zone: while the 
planet might be unsuited to the evolution of life, its moons 
might be. Moons with masses of at least one third Mq, 
and orbiting around gas giant planets in the habitable zone 



may have habitable environments (Williams et al. 19971 



The moon's orbital stability plays a role in habitability as 
well. Clearly, if the moon leaves orbit, it will probably 
leave the habitable zone. 

Although extrasolar moons have not yet been found, 
several methods to detect them have been investigated. 



After Kaltenegger ( 2010 1 , the following methods can detect 



extrasolar moons: 



1. Transit timing variations (Sartoretti & Schneider 
T999] [Agol et al.(2005) |Holman fc Murray||2005f ~ 



2. Transit duration variations (Kipping 2009). 

3. Light curve distortions ( Szabo et al.|2006 ). 

4. Planet- moon eclipses ( Cabrera fc Schneider|2007 ). 

5. Microlensing ( |Han||2008 1 . 

6. Pulsar timing ( Lewis et al.p 008). 



7. Distortion of the Rossiter-McLaughlin effect of a 
transiting planet ( Simon et al.||2009 ). 



Considering the speed at which observational instrumen- 
tation has developed, it is only a matter of time before 
extrasolar moons are discovered. 

Tidal torque is important to the long-term orbital sta- 
bility of extrasolar moons. A binary system can be in 
tidal equilibrium only if coplanarity (the equatorial planes 
of the planet and moon coincide with the orbital plane), 
circularity (of the orbit), and corotation (the rotation pe- 
riods of the planet and moon are equal to the revolution 
period) have been fulfilled. Further, stability occurs only 
if the orbital angular momentum exceeds the sum of the 
spin angular momenta of the planet and moon by more 
than a factor of three flHut||1980| ). 



Counsclman (1973) studied the stability of these equi- 



libria only with respect to coplanarity and circularity. He 
pointed out that in a planet-moon system with lunaij^] 
tides, there are three possible evolutionary states. 

Counselman state (i):the semi- major axis of the moon's 
orbit tidally evolves inward until the moon hits the planet. 
Example :Phobos around Mars. 

Counsclman state (ii):the semi-major axis of the moon's 
orbit tidally evolves outward until the moon escapes from 
the planet. No solar system examples are available. But 
this result would happen to the Earth-Moon if Earth's 
present rotation rate were doubled. 

Counselman state (iii)dunar orbital and planetary spin 
angular velocities enter mutual resonance and are kept 



commensurate by tidal forces. Example:Pluto and Charon. 
This Counselman state is static, while state (i) and (ii) are 
evolutionary. 

Here, we consider a star-planet-moon system with stel- 
lar tides. Although they did not consider the effects of 
lunar tides or maximum distance from the planet, Ward & 



Reid ( 1973 ) examined the impact of solar tides on plane- 



tary rotation in a limited star-planet-moon system. 



Barnes & O'Brien (2002) considered a similar case, con- 



sidering the maximum distance of the moon but neglecting 
the lunar tide's effect on planetary rotation. They found 
just two possible final states:the moon may either hit the 
planet or escape from it. 

In this paper, we consider a star-planet-moon system 
with both stellar and lunar tides, and lunar maximum dis- 
tance from the planet. Stellar and lunar tides both affect 
planetary spin, whereas stellar and lunar tides affect planet 
and moon orbits, respectively. We do not consider the ef- 
fect of stellar tides on the moon's rotation. Stellar tides 
should sap angular momentum from the system but this 
effect is less important if the mass of the planet is at least 
ten times greater than the mass of the moon. We apply 
tidal theory and set up a system of differential equations 
that govern the planetary rotational rate and orbital mean 
motion as well as the orbital mean motion of the moon. 
The system of differential equations is solved numerically. 
Finally, a formula for the length of time the moon will 
stably orbit is found. We then apply this result to hypo- 
thetical extrasolar planet moon systems. 

2. THEORY 

In this paper, we use standard tidal evolution theory 



with the constant Q approach (Goldreich & Soter 19661, 



3 In this paper, we use "lunar" as the adjective of any moons, not just 
the Moon 



along with the following assumptions: 

1. The spin angular momentum of the planet is parallel 
to the orbital angular momenta of both the moon 
about the planet and the planet about the star; i.e. 
the planet has 0° obliquity, the moon orbits in the 
planet's equatorial plane, and the planet and moon 
motions are prograde. 

2. The total angular momentum, that is the sum of the 
moon's orbital angular momentum and the planet's 
rotational and orbital angular momenta, is constant. 
We neglect the orbital angular momentum of the 
moon about the star and the moon's rotational an- 
gular momentum. 

3. The moon's orbit about the planet and the planet's 
orbit about the star are circular. 

4. The moon is less (at most ~ 1/10) massive than the 
planet and the planet is also less (at most ~ 1/10) 
massive than the star. 

5. The star's spin angular momentum is not considered 
nor are the planet's tides on the star or the star's 
tides on the moon. 
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6. The specific dissipation function of the planet, Q p , is 
independent of the tidal forcing frequency and does 
not change as a function of time. 

Planetary Q p falls into two groups. The first group 
has values of Q p that range from 10 to 500. The ter- 
restrial planets and satellites of the Jovian planets are in 
this group. The other group has Q p values greater than 
10 4 . The Jovian planets are in this group. In the case of 
the Earth, tidal dissipation is due to friction between the 
tidally generated currents and the ocean floor and occurs 
mostly in shallow seas. For Mercury and Venus, tidal dis- 
sipation is driven by viscous dissipation within the bulk 
planetary interior. The mechanism for tidal dissipation 
within giant planets remains unknown. 

In that case, the torque on the planet due to the moon 



Tp- m is given by 


Barnes & O'Brien 


(2002 


) Goldreich & 


Soter 


(1966 


); 


Murray & Dermott ( 


2000 1 in Chapter 4 



' p—rn 



3 k 2p GMlRl 

-sgn(S2 p 



?P a tn 



(1) 



where tt p is the rotational rate of the planet, k 2p is the 
tidal Love number of the planet, G is the gravitational 
constant, R p is the radius of the planet, M m is the mass 
of the moon, a m is the semimajor axis of the moon's orbit, 
and n m is the orbital mean motion of the moon. The 
function sgn(J7 p — n m ) is 1 if (il p — n rn ) is positive, -1 if 
(£lp — n m ) is negative, and undefined if (Jl p — n m ) = 0. 

Similarly, the torque on the planet due to the star t p ^ s 



is 



3 k 2p GM*R 5 p 



n p ), 



(2) 



where M s is the mass of the star, a p is the semimajor axis 
of the planet's orbit, n p is the orbital mean motion of the 
planet. 



For the spin angular momentum of the planet, 



dfl p 



dt 



' p—m 



' p—Sf 



(3) 



where the planet's rotational moment of inertia I p = 
aM p R p . a is the moment of inertia constant. For Earth 



and Jupiter, ct's are 0.3308 and 0.254, respectively (de Pa- 
ter fc Lissauer|2001 |. 



The change in orbital angular momenta of the moon 
about the planet and the planet about the star, by New- 
ton's Third Law, are equal and opposite the moon's and 
star's torques on the planet, respectively: 



and 



where L„ 



dL m 
~dt~ 

dL p 
~dt 

,2 



Tp—s j 



(4) 



(•5) 



= M m a l m n m and L s = M p a p n p . 



Using Kepler's Third Law, n^a^ w GM p and n p a p » 
GM S because we assume that M m <C M p and M p <C M s . 
This allows us to derive these expressions for L m and L p 



M m (GM p f/ 3 



1/3 



and 



M p (GM s f/ 3 



1/3 



(6) 
(7) 



To determine how n rn and n p vary with time, we take 
the derivative of L m and L p with respect to t, set the 
results equal to equation Q and equation ([5]), and solve 



for ^ ^d ^ 
Then, we have 



dn r , 



3r ; 



and 



dt 



drip 
~dt 



p-m n 4/3 



M m {GMp)Vz 



3t„ 



,4/3 



M p {GM s )Vz p 



(8) 



(9) 



When we combine equation ([T]), equation ([2]), equation 
([3]), equation ([8]), and equation ([9]), we obtain the dif- 
ferential equations that govern the time-evolution of the 
star-planet-moon system: 



< 



dn m _ 9 k 2pK GM„ 16/3 /p N 

dn p _ 9 k 2p Rl 1 r7 16/3 SffnCO -71 \ 

dt — 2 Q p (GA/p)(GM 3 ) 2 / 3 "V b & 1 H"p ll P> 

dn p _ 3 fc2 P flp (GM m f 4 f p _ v 

dt ~ 2 Qp a(GM p ) 3 ''m°& u V si P "m) 

~2 Q P P a(GMp) n p S S n (^P _ n p)- 



(10) 



The solutions to these differential equations are 



n m (t) 



n p (t) = 



n p (t) 



39 k 2v Rl GMr, 



2 Q p (GM p )8/3 



t sgn(O p - n m ) 



+n- 13 / 3 (t = 0))(lli) 
t sgn(i7 p — n p 



39 k 2p R 5 p 



1 



! Qp (GM p )(GM s ) 2 /3 
1 f GM, 



-3/13 



+V 13 / 3 (o)) ' (lib) 



+O p (0). (11c) 



These solutions are only valid if each of sgn(fi p — n m ) 
and sgn(f2 p — n p ) are constant in time. Also, these so- 
lutions are only valid when the planet's rotation is not 
tidally synchronous with either the star or the moon, i.e. 
il p — n m and f2 p — n p ^ 0. When the planet's rotation 
is synchronized, we must use an alternative approach. 
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When synchronization has occurred, i.e. when £l p = n m 
or fl p = n p , we follow the evolution of the system using 
conservation of angular momentum: 



M m (GM p ) 2 / 3 



aR 2 M p n p (t)- 



M p (GMt 



\2/3 



= L , (12) 



where L = ^ /3( gJ + aR z p M p Vl p {Q) + ^ /3( ^ is 
the initial angular momentum. 

By our assumption 2, the total angular momentum is 
the sum of the moon's orbital angular momentum, which is 
the first term, the planet's rotational angular momentum, 
which is the second term, and orbital angular momentum, 
which is the third term. 

When the planet is not tidally locked with either the 
star or the moon, these three equations are valid: 



,(*) = 



n p (t) = 



39 k 2p R 5 p GM m 
T Q p (GM p )8/3 



t sgn(f2 p 



+n m 13/3 (0)r 3/13 (13a) 



1 t8m(n -n) 

2 Q p (GM P )(GM 8 )V* m p p) 

\ -3/13 

+v 13/3 (°)) ( 13b ) 



(GM m )(GM p ) 2 / 3 



+ 



aR 2 p (GM p )fl p (t) 

(GM P )(GM S ) 2 / 3 
n p /3 (t) 



GL a . (13c) 



Even though equation (11c) and equation (13c 



are 



equivalent, equation (13c) is valid when the planet's rota- 



tion is tidally locked to the moon because equation (13c 



is derived from the conservation of angular momentum. 
When the planet is tidally locked with the moon, i.e. 



= fip, equation (13a) is not valid. Hence, in that case, 



n p (t) 



'39 k 2p R 5 p 



1 



2 Qp~ (GMp)(GM s )V3 sgn ^ m " n P> 

-3/13 



+n' 



-13/3 



(0) 



(14a) 



(GM m )(GM p ) 2 / 3 
n m (t) 



+ aR 2 p {GM p )n m (t) 

| (GM P )(GM S ) 2 / 3 
n p /3 (t) 



= GL . (14b) 



When the planet is tidally locked with the star, i.e. n p 



f2p, equation (13b) is not valid. Hence, in that case, 



n m (t) 



'39 k-2 P Rl GM m 
~2 Q p (GM P )8/ 3 



t sgn(n p - n m ) 



-13/3 



(0) 



-3/13 



(15a) 



(GM m )(GM p ) 

Urk (t) 



2/3 



aR 2 p {GM p )n p {t) 



(GM P ){GM S ) 
n p /3 (t) 



GL Q . (15b) 



3. Numerical Solutions 



We first explore the implications of equations ( 13 ), ( 14 ), 



and ( 15 ) numerically. In these simulations, we start with 
f2p(0) > n m (Q) > n p (0). Physically, this condition implies 
that one planet year is longer than one planet day and that 
the orbital period of the moon is between them. A typical 
example is our Sun-Earth-Moon system. The shapes of 
the resulting graphs of ft p , n m , and n p as a function of 
time depend on the torques due to the planet and moon, 
but also on the orbital angular velocities of the planet and 
moon. If both the orbital angular velocities of the moon 
and planet are slower than the spin angular velocity of 
the planet, then both the torques due to the moon and 
those due to the star brake the rotation of the planet. If 
the orbital angular velocity of the moon is faster than the 
spin angular velocity of the planet, then the spin angular 
velocity of the planet may increase or decrease, depending 
on the relative magnitude of the torque due to the moon 
and star. 

One important lunar escape condition for the calcula- 
tion is the critical semimajor axis; this is the outermost 



stable orbit for the moon. Barnes & O'Brien (2002) stated 



that the critical semimajor axis, a cr u, is 

O-crit = fR-H, 



(16) 



where / is a constant and Rr is the radius of the Hill's 
sphere ( de Pater k. Lissauer||2001 ) 



R 



H 



1/3 



(17) 



where a p is the semimajor axis of the planet. Orbits out- 
side a cr it are not stable. In this paper, we follow the orbits 
using angular velocity instead of semimajor axis. Consid- 
ering n cr it such that 



l crit a crit ~ G-Mp, 



(18) 



then, using (16) and (17), and Kepler's Third Law for the 



planet, n p a p — GM S , we calculate n cr n as a function of /: 



n crit{t) — 



1/2 



^p(^). 



(19) 



The value of / is not well-determined. Barnes fc O'Brien] 



(2002) used / = 0.36. Domingos et al. (20061 suggested 



/ = 0.49. In this paper, we use / = 0.36 for numerical 
calculations because it is the most conservative estimate 
for the moon to remain bound. 

This is the critical mean motion that is the lowest stable 
angular velocity for the moon. From section [4T] to section 
4.31 we enforce that 



(t) < n m (t). 



(20) 
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This means that the moon has a stable orbit. 

In the resulting numerical integrations, we found three 
classes of stable outcomes. We call the three stable out- 
comes : 



• Type I(Fig[T]and[2} 

— planet-moon become synchronous 

• Type II(Fig§ 

— planet-star become synchronous first, then 
planet-moon become synchronous later 

• Type III(Figj4]) 

— planet-moon never synchronous. 

For Type I, there are two subcases. 

In each of the three stable outcomes, the first part is 
common. Initially, f2 p (0) > n m (0) > n p (0). Since the 
orbital angular velocities of the moon and the planet are 
slower than the spin angular velocity of the planet, the 
torques due to the moon and the star brake the rotation of 
the planet. This continues until the spin angular velocity 
of the planet is equal to the angular velocity of the moon. 
We call this time Tl. From the beginning to Tl, the planet 
loses rotational angular momentum to the orbital motions 
of the moon and the planet. By gaining angular momen- 
tum, the orbital motions of the moon and the planet slow 
down and their semimajor axes increase. 

After Tl, each type has its own characteristics. An- 
other feature that Type I, II, and III have in common is 
that the planet's angular velocity, n p (t), always decreases 
due to solar tides. This indicates that the orbital angular 
momentum of the planet always increases. 

In Type I, if the tidal locking starts at Tl, then the 
system is Type I Case 1 (FigQ. The system is Type I 
Case 2 if the tidal locking starts after Tl (Fig{2|. 



Type I 



Casel (Fig[T]) 

* In our Type I Case 1 star-planet-moon sys- 
tem, the torque on the planet due to the 
moon is greater than that due to the star 
at t = Tl. At Tl, the planet and the 
moon then assume a synchronized state 
with Q p = n m . Once they reach this syn- 
chronized state, they will stay in this state 
until the end for Type I Case 1. Since 
the tidal torque on the planet due to the 
moon is greater than that due to the star, 
the moon's orbital velocity, n m , and the 
planet's spin angular velocity, fl p , are kept 
equal. In this synchronized state, only the 
orbital motion of the moon loses angular 
momentum; the planet's orbital and spin 
motion gain angular momentum. 



«m(t) 
Cl p (t) 
n p (t) 



3x10'° 4x10" 
Time (years) 



5x 10 10 



Fig. 1. — Here we graph the time evolution of £l p , n p , 
and n m for Type I Case 1. We use the present data 
of our Sun-Earth-Moon system for the initial condition, 
i.e. n m (0)=84 rad/year, fl p (0) — 730n rad/year, and 
n p (0) = 2ir rad/year. We ta ke ki v and Q p for E arth to 
be 0.299 and 12 respectively ( |Murray fc Dermott| ( |2000| ), 
pgl66). The black vertical line on the right corresponds to 
the lunar orbital frequency at which the moon is orbiting 
at the planetary radius, i.e. when it crashes into the Earth 
and is destroyed. 

- Case2(Fig[|) 

* In Type I Case 2, the moon's tidal torque 
on the planet is slightly smaller than the 
stellar torque at Tl, but the planet's ro- 
tation never becomes tidally locked to the 
star. There is a brief period when tt p is 
between n m and n p , tidally locked to nei- 
ther the star nor the planet. At t = Tl, 
the planet and the moon cannot reach the 
synchronized state because the torque due 
to the moon is smaller than that due to 
the star. The planet's spin keeps decreasing 
for a while. In this period, the moon's or- 
bital motion and the planet's spin motions 
lose angular momentum, and the planet's 
orbital angular momentum increases be- 
cause of the decreasing semimajor axis of 
the moon. As the moon's orbital angular 
velocity, n m (t), increases, so does the tidal 
torque due to the moon. Shortly thereafter, 
the torque due to the moon overcomes the 
torque due to the star. The planet's spin 
angular velocity, ft p (t), starts to increase. 
Then, the planet and the moon reach the 
synchronized state. Once synchronous, the 
moon's orbital motion loses angular mo- 
mentum, and the planet's orbital and spin 
motions increase angular momentum. This 
case is distinct from Case 1 in that there is a 
period when the moon is migrating inward, 
but is not synchronized with the planet's 
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Time (years) 

Fig. 2. — This graph is for Type I Case 2. We use the 
same conditions as in Fig [l] except the planet's mass is 1.2 
M®. Note the "notch" in fi p just after Tl; it is what 
differentiate Case2 from Casel. 



spin. 

Type II(Fig§ 



and the planet's spin motion both lose angu- 
lar momentum. At the star-planet synchronized 
state, only the moon's orbital motion loses an- 
gular momentum. 

We can sec the difference between Type II and 
Type III in Stage 3. Roughly speaking, if we can 
see Stage 3, then the system is Type II. If Stage 
3 is so short that we cannot see it, the system 
is Type III. For Type II, the torque due to the 
star becomes smaller than the torque due to the 
moon as the moon spirals inward. The planet's 
rotation becomes tidally locked to the moon, 
after which only the moon's orbital motion loses 
angular momentum. 

• Type III(Figj4]) 

— In this case, the tidal torque due to the star 
is always greater than that due to the moon. 
The amount of loss or gain in angular momen- 
tum for the moon's orbital motion is so small 
that we can treat the sum of the orbital an- 
gular momentum and the spin angular momen- 
tum of the planet as a constant. In essence, the 
planet's spin evolves as if the moon does not ex- 
ist - this corresponds to the |Barnes &: Q'Brien| 
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Fig. 3. — This graph is for Type II. We use the same con- 
ditions as in Fig{T] except the planet's mass is 4 M®. Here 
the planetary spin becomes synchronous with its orbit af- 
ter Tl for a time T2. But spins up to become synchronous 
with the moon thereafter. 



lxlO 10 2xl0 10 3xl0 10 4xl0 10 5xl0 10 

Time (years) 



Fig. 4. — This graph represents Type III. We use the same 
conditions as in Fig[I] except that the moon's mass is set 
to 0.1 Mj. The planetary rotation becomes synchronous 
with the star, and never with its moon. 



— For Type II, at t — Tl, the tidal torque due to 
the star is greater than that due to the moon, 
which forces the planet's rotation to continue to 
slow down until it becomes synchronized to the 
star (Tip = rip). 

In Stage 2, then, the planet and star remain 
in a synchronized state because the torque due 
to the moon does not overcome the torque due 
to the star. Until the planet and star reach a 
synchronized state, the moon's orbital motion 



All n p s in Fig{TJ-Fig|4]decrease only a very small amount. 
These changes are almost unnoticeable. However, the dif- 
ferences between the outcomes of tidal evolution in a two 
body system and a three body system come from these 
small changes. Because n p decreases in our model, which 
does not include the star's tidal response to the planet, 
our results have a systematic error with respect to the life- 
time of moons for close-in planets which experience orbital 
decay. 
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4. Analytical Lifetimes 

In this section, we derive analytical formulae for the 
total moon lifetime for each of Type I, II, and III. The 
total moon lifetime means the time it takes for the moon to 
either hit the planet or escape. We will discuss the results 
here; The full derivations of the formulae are described in 
the Appendix. 



escape from the planet. At the synchronized state, the 
orbital velocity of the moon is equal to the spin angular 
velocity of the planet, i.e. n m — tt p . This corresponds 
to Counselman state (hi). After the synchronized state, 
the moon has a brief spiral inward orbit (Fig [5| ; this is 
Counselman state (i). 

4.2. Type II Solution 



4.1. Type I Solution 

As discussed in section 3, Type I has two different cases. 
Because we can use the same formula to calculate the life- 
time of the moon in each case though, we call both cases 
Type I. By creating a new function, n m (t), which coincided 
with n m (t) after Tl, we found the formula for the lifetime 
of the moons. 



The formula for the lifetime of the moons for Type I, T, 



is: 



Qp 



:(GM P )(GM S ) 2 / 3 
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3 3 / 4 (GM p )(GM 3 ) 2 / 3 
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After the synchronized state is broken at the very end 
of the moon's life, the moon has a spiral inward orbit (Fig 
[5]). We do not explicitly include this period because it is 
very small compared to T. For Type II, we have the same 
situation. 



> 



For Type II, there are three stages (Figj3j). We calcu- 
lated the time intervals for Tl, T2, and T3, respectively. 
By adding them up, we found the lifetime of the moons for 
Type II. 

The formula for the lifetime of the moons for Type II, 
T, is: 



T = Tl + T2 + T3 
— 2T1 1 2 



39 k 2p R p 
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(3 3 / 4 GL -4{(GAf m ) 3 (GM p ) 3 a^} 1/4 ) 3 

3 39 / 4 (GM p ) 12 (GM 3 ) 8 
(GLp) 13 



(22) 



{(GM p y/ 2 (GM m y/ 6 + (GM p )(GM s ) 2 / 3 } 
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Fig. 5. — This graph is a magnification of the last part of 
Fig[l] At the very end, the synchronized state is over. The 
orbit of the moon decays inward much faster than before. 
Because the duration of this final death spiral is so short 
- 7000 years in a 67.5 Gyr evolution - we neglect it in our 
analytical formulations. 

Interestingly, a Type I star-planet-moon system experi- 
ences all three Counselman states. From t = to Tl, the 
orbital velocity of the moon, n m , decreases. This indicates 
that the orbital semimajor axis of the moon increases. This 
is Counselman state (ii), except that the moon does not 



We could not find the analytical expression for Tl. 
However, we can calculate Tl by solving the following sys- 
tems of equations numerically for t 



(t)= (39 ^pK GM m t , -13/3/ Q ^ 
«mW ~ \ 2 Q p (GAfp) 8 / 3 t + Um v u J 



-3/13 



n\ = ( 3g k 2pR'p i 

'WJ \ 2 Q p (GM„)(CM,)'/3 



M ^ 2/3 +aHf ) M p n m (t) 



t + n- 13/3 (0)) 

M P (GM 3 ) 2/3 __ 



-3/13 



(23) 



n P /3 (t) 



Tn 



After Stage 3, there is a brief Stage 4 wherein the moon 
makes its final death spiral into the planet's cloud tops. 
At Stage 4, n m (t) ^ O p (t) - actually, n m (t) > O p (t). Since 
T4 is very small compared to Tl, T2, and T3, we do not 
explicitly include T4 in our calculation. 

Type II has all three Counselman states, like Type I, 
plus one extra state. Stage 1 corresponds to Counsel- 
man state (ii) except that the moon does not escape from 
the planet. Stage 2 is the extra state. At this stage, the 
planet and the star are tidally locked. Because Counscl- 
( 1973 ) considered a planet-satellite system, the planet 



man 



could not be tidally locked with the star, hence Stage 2 has 
no corresponding Counselman state. Stage 3 corresponds 
to Counselman state (iii) . The planet and the moon reach 
a synchronized state. At Stage 4, the moon has a spiral 
inward orbit. This is Counselman state (i). 

4.3. Type III Solution 

For Type III, we can calculate the lifetime of the moon 
using symmetry (Fig 19) as for (Barnes & O'Brien 2002). 



The formula for the lifetime of the moons for Type III, 
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T, 



is: 



T = 2T1 



2 Q p {GM l 



N8/3 



39 k 



2p 



Rl 



GMr, 



n m 13 / 3 (0). (24) 



In general, the system is Type III if the moon is very 
small compared to the planet. If we set (GM m ) 3 {GM p ) 3 aR p 



and (GMp^^tGMm) 7 / 6 in equation (22) equal to zero, 



then we get equation ( 24 ) 



Type III has one Counselman state and one extra state. 
The first part, from t = to Tl, is Counselman state (ii) 
except the moon does not escape from the planet. The 
second part, t > Tl, is the extra state which is the same 
as Stage 2 in Type II. 
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4.4. Type IV Solution 

So far, we assume that the orbit of the moon is always 
stable. In this section, we calculate when the orbit of the 
moon becomes unstable. This means 



(25) 



n crit (t) > n m (t) 



at some t (Fig|6j). In this case, we can use equation ( pL3| ) 
with sgn(f2 p — n m )—l and sgn(f2 p — n p )=l because the 
planet is not tidally locked with either the star or the 
moon. To find the time when the orbit of the moon be- 
comes unstable, we set n cr i t (t) = n m (t). Then we solve for 
t. The lifetime of the moon in this case is 



T 



2 Q P ( (GA/ P ) 8/3 (GM S ) 2/3 

39 k 2v Rl \ (GM m )(GM s )2/3-(/3/ 3 )i3/6(GAf p )V3 , 

x ((/ 3 /3) 13 /V 13/3 (0) - n„ 13/3 (0)) . 



Type IV has just one Counselman state: Counselman's 
state (ii). 

5. Determine the Type of the System 

The expressions in section 4 can be used to calculate the 
ultimate lifetime of any star-planet-moon system, provided 
you know which type of system it is. In this section, we 
show how to determine a system's type. 

5.1. Condition for Type I Case 1 

The condition for Type I Case 1 is that the magnitude of 
the torque due to the moon is greater than the magnitude 
of the torque due to the star at t = Tl (FigjTJ) , 



\r P - m (Tl)\ > |r p _.(Tl)|. 



(27) 



This condition implies 



Tl < 



2 Q p (GM p )(GM m ) 7 / 6 (GM S ) 2 / 3 
39 k 2p Rl (GM P )V 2 (GM S )V 3 - (GM m )V6 

( GM P ^ 13/6 



x < n 



-13/3 
P 



(0) 



V GM, 



£) n - 13/3(o) • (28) 



Fig. 6. — This is the graph of n m (t) and n cr i t (t) for Type 
IV. The dashed line is n cr it{t) and the solid line is n m (t). 
At t = T, the moon is on the outermost stable orbit. After 
that, the orbit becomes unstable. Then, the planet loses 
the moon to interplanetary space. In this case, the life- 
time of the moon is T. We use the data of the present 
Sun-Earth-Moon system except with an initial Earth spin 
angular velocity, O p (0), of 12007T. 

This is the condition for the Type I Case 1. We can get 



Tl by solving equation ( 23 ) numerically 



If the system satisfies equation (28 1, then we can con 



elude that it is of Type I Case 1. If not, the system may 
be Type I Case 2, Type II, or Type III. 



The sign on the right side of equation (|28|) depends on 

(29) 



-13/3 



(0) 



/ GM P ^ 13/6 



\GM„ 



-13/3 



(0) 



because (GM P ) 1 / 2 (GM S ) 2 ^ 3 — (GM TO ) 7 / 6 > (10 11 / 6 



l)(GM m ) 7 / 6 > by our assumption 4. If equation (29) 
is negative, then the system cannot satisfy equation ( 28 1 . 



Hence, the system is Type I Case 2, Type II, or Type III. 



The inequality 



-13/3 



MS) 



13/6 



n™ 3/3 (0)!> : tl 



implies that |r p _ m (0)| < |r p _ s (0)|. This means that if 
the initial torque due to the star is greater than the initial 
torque due to the moon, the system cannot be Type I Case 
1. 

5.2. Condition for Type I Case 2 

Assume n p (Tl) and n TO (Tl) are known. Let t + be the 
time from Tl, when the magnitudes of the two torques are 
equal (Fig[7|). 

The condition for Type I Case 2 is that 

fl p (t + ) > n p {t + ) (30) 

where t t satisfies 

\r p - m (t+)\ - |r p _ a (i+)|. (31) 
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Fig. 7. — This graph is a magnification of Fig [2] The 
rotational rate of the planet, f2 p (i), decreases until t = 
Tl + 1 + because the torque due to the star is greater than 
the torque due to the moon. At t = Tl + t + , these two 
torques are equal. After that, the torque due to the moon 
exceeds the torque due to the star. The rotational rate of 
the planet, Q p (t), starts to increase. Then, the planet and 
the moon reach a synchronized state. 



Fig. 8. — This is a graph of Type III after the synchronized 
state of the star and planet end. In this graph, t = 
means the time when the synchronized state of the star and 
planet ends. We used the present data of our Sun-Earth- 
Moon system except M p = 18M®, n m (0) = 240, n p (0) = 
fip(0) = 6.28 for n m {t). For n m (t), we used n p (0) = 6.28, 
and n m (0) = O p (0) = 343. 



This condition implies that 

ai b 12 X 4 - GL X 3 + % 3 < 0, 



where 

r — - 1 

a R p (GM p ) 

ai = (GM^/\GM m f/ 7S 

OL2 = 



(32) 
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b = {(GM m ) 7/6 + (GM P ) 1/2 (GM a ) 2 '' 3 } 

{(S) 13/6 ^ 13/3 (") + (GM g^ )a/ ' np- 13/3 (ri)} 



X 



If the system is not Type I and satisfies equation ( 32 1 , then 
it is Type I Case 2. 

5.3. Conditions for Type II and III 

As you can see in FigjSJ the planet and the moon reach 
a synchronized state at Stage 3 for Type II. In Fig(3J the 
synchronized state of the star and planet seems to continue 
indefinitely. The planet and star seem to be tidally locked 
until the end. Actually, it ends when the moon spirals into 
the planet, which occurs when n m (t) is large enough. In 
Fig|4j the total lifetime of the moon, T, is about 53.3 billion 
years. The synchronized state of the star and planet ends 
at about 52.7 billion years. In Type III, there is a stage 
that corresponds to Stage 3 for Type II. However, because 
that stage is short compared to the total lifetime, it is hard 
for us to see it. 

Fig(H]is a graph of Type III after the synchronized state 
of the star and planet ends. The planet is so big that 
the spin angular velocity of the planet cannot increase fast 



1/13 



enough to catch up to the orbital angular velocity of the 
moon. n m (t) is the hypothetical situation in which the 
planet and the moon are tidally locked from the beginning. 
We introduced n m {t) in 



4.1 



Fig|9] is a graph of Type II at Stage 3. The planet 
is not big enough so that the spin angular velocity of 
the planet increases fast enough to catch up the orbital 
angular velocity of the planet. n m (t) is the graph of 



39 k 2p R% QMn 



n- U/3 (Tl 



\ 2 Q p (GM p )s/3 

same equation of n m (t) in Type III. 



-3/13 



that is the 



T5 is the maximum range of n m (t) and T6 is the max- 
imum range of n m (t) in Type II and n m (t) in Type III. 

From Fig|8]and Figj9) the condition for Type II is T5 > 
T6 and the condition for Type III is T5 < T6. 

The condition for Type II, T5 > T6, implies that 



(3 3 / 4 GL - 4 {(GM m f{GM p faRl) llA 
((GM m ) 7 /6 + (GAf p ) 1 /2( GMs )2/3 
> 3 39 / 4 (GM p ) 6 (GM s ) 8 (GL ) 13 . 



V ) 
, 12 



13 



(33) 



Similarly, the condition for Type III, T5 < T6, implies 
that 



(3 3 / 4 GL - 4 {{GM m f(GM p faRl) 1/A 
((GM m ) 7 / 6 + (GM p )V2( G M 5 .)2/3) 12 
< 3 3g / 4 (GM p ) 6 (GM s ) s (GL ) 13 . 



13 
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and usually much smaller. The fractional error is defined 
as the absolute value of the lifetime from simulation minus 
lifetime from formula divided by lifetime from simulation. 

In the first example, Q p =12 is not realistic for a high 
mass planet. However, we want to show just how the mass 
of the planet affects the lifetime of the moon. To get a 
more realistic result, we should change both ki p and Q p 
for Jovian planets. If we did that, we would expect that 
the lifetime of the moons would be much longer. 

In both examples, the lifetime in most of these figures 
is longer than the typical main sequence lifetimes of solar- 
type stars. Over such a timescale additional processes, 
such as inflation of the stellar radius, and the resulting 
changes to tidal evolution become important. 



Fig. 9. — This is a graph of Type II at Stage 3. In 
this graph, t = means t = Tl + T2 in the original 
Type II graph. We used the present data of our Sun- 
Earth- Moon system except M p = 15M®, n m (0) = 219, 
n p (0) = fi p (0) = 6.28. For n m (t), we used n p (0) = 6.28, 
and n m (0) = Q p (0) = 286. 

5.4. Condition for Type IV 

Looking at the graphs of Type I (FigjTJ Figj2]), Type II 
(Fig|3]), and Type III (Figg, we know n. m (t) has a mini- 
mum at t = Tl. Hence, the condition for Type IV is 



Tl > T 



Qp 



(GM p f 3 (GM s f 3 



39 k 2p Rl \^(GAf m )(GM a ) 2 /3_(/3/ 3 )i3/6( GMp )5/3 



((/ 3 /3) 
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n m 13/3 (0) 



(35) 



We can find Tl by solving the system of equations ( 23 ) 
numerically. 



6.1. 



First Example - Changing The Mass of The 
Planet 



Fig 10 shows how the lifetime of a 1 Mj moon varies as 
the mass of the planet changes from 0.1 to 2 M ffi . Fig 



11 



similar, but with the mass of the planet varying from 2 to 
25 M®. The difference between the data from numerical 
simulations and the graph generated by the formula is very 
small. From the graph FigfToJ and FigrTT] every data point 
is on the curve. The transition from one type to another 
is smooth. 



As you can see in Fig 10 and Fig 11 the lifetime of the 
moon increases as the planet mass increases. The increased 
longevity occurs because the effects of the lunar and stellar 
tides on the planetary spin evolution are reduced as we 
increase M p , i.e. the planet is not braked as easily by each 
of the these effects as we go to more massive planets. The 



system therefore continues to evolve with 51 p > 



for 



longer, in which the tidal torque on the moon is positive, 
lengthening the overall evolution. 



6. Applications 

In this section, we check the formulae for the lifetime 
of the moon and the condition for the type of system first 
in two examples, and then showing an application of our 
results. 



In the first two examples, we use the present data of our 
Sun-Earth-Moon system. We take k2 P and Q p for Earth to 
be 0.299 and 12, respectively ( |Murray fc Dermott| (2000| ), 
pgl66). In the first example, we survey a range for the 



6.2. Second Example 
Moon 



Changing The Mass of the 



mass of the planet from 0.1 to 25 M®. In the second ex- 
ample, we explore moon masses from 0.01 to 2 Mj. For 
Type I, we can plot the graph of the lifetime of the moon 
easily because every parameter is constant (Eq|2l"|). For 
Type II and III, we need to know the expression for Tl to 
plot the graph (Eq 22|24 ). In order to have Tl, we must 
calculate Tl numerically. Then, by using an approxima- 
tion method, such as best-fitting or interpolation, we can 
obtain the expression of Tl. 

The fractional error of our analytical formulae is always 
smaller than 10~ 3 when compared to numerical integration 



Fig[l2] shows how a moon's lifetime would change as 
the moon's mass varies from 0.01 to 2 Mj. For low-mass 
moons, the greater the mass of the moon, the more quickly 
the moon evolves from equation (10). This is why the 
heavier moon has the shorter lifetime in Type III (Barnes 
fc O'Brien|2002| ). Additionally, at Type I Stage 1, a heavy 
moon evolves faster than a light moon. Indeed, Tl de- 
creases as the mass of the moon increases. 

However, once the moon and the planet reach their syn- 
chronized state, the planet keeps the heavier moon for a 
longer time. When the planet and the moon become tidally 
locked, this planet-moon system behaves like one object. 
As the mass of the moon increases, the moment of iner- 
tia of the planet-moon system increases. Because the star 
saps angular momentum at a constant rate, the planet- 
moon system evolves more slowly when the system has the 
heavier moon. 

Between Types III and I, the lifetime of the moon has 



10 




0.5 1.0 1.5 

Mass of the planet (M ffi ) 



5s. 

1 1x10 



13 



§ 8xl0 12 
o 

« 6xl0 12 
-g 

I 2xl0 12 





5 10 15 20 
Mass of the planet (M ffi ) 



25 



Fig. 10. — This graph shows the lifetime of a hypotheti- 
cal system with Sun-Earth-Moon parameters with varying 
planetary mass from 0.1 to 2 AT®. Each dot represents the 
result from numerical solutions, and the curve is generated 



by equation (21 1 or equation (22), depending on the type 



of system. The lifetime of the moon increases linearly as 
the mass of the planet increases in this region. The bor- 
ders between Case 1 and Case 2, and Case 2 and Type II 
are M p = 1.04M® and M p = 1.27M®, respectively. The 
border between different types depends on the initial con- 
ditions and Q p . 



a minimum value. When the mass of the moon is 0.307 
Mj, it has a minimum lifetime of 4.24 x 10 10 years. 

We know that the tidal effect of the real Moon is impor- 
tant on the Earth. To show the utility of our approach, we 
calculate the lifetime of hypothetical moons with and with- 



out lunar tidal effect (Fig 13 1. For low masses of the moon 



there is no difference between these results and previous 



results of Barnes & O'Brien (20021 because the effect of 



the lunar tides is small. For the high mass moons, the life- 
time of the moon with the lunar tidal effect is significantly 
longer than that without lunar tidal effect. 

6.3. Third Example - Other Systems 

Finally, we apply our results to extrasolar star-planet- 
moon systems where we expect that the first exomoons will 
be discovered. To see the big picture, we chose some typical 



combinations of the stars, planets, and moons (Fig 14 ). We 



choose 1.0 Mq and 0.3 Mq stars as the parent stars. For 
each parent star, we investigate 7 planet/moon systems. 
For rocky planets, we chose Earth/Moon, Earth/Mars, 
and 8-A/0-super-earth-planet/Venus systems. For ice gi- 
ant planets, we choose hypothetical Neptune/Earth sys- 
tem. For gas giant plants, we choose Saturn/Earth, 
Jupiter/Earth, and 10-Mj up -planet/Earth systems. To see 
the trend of each type of the planets, we also separate the 
systems by Q p (FigjlH). We use Q p = 100, 10 4 , and 10 5 
for rocky, ice giant, and gas giant planets, respectively. 

Fig [14] shows the moon stability lines for 1 to 10 Gyr ap- 
plied to types of planet/moon systems. It is worth noting 
that a star of mass 0.3 Mq has a main sequence lifetime 



Fig. 11. — Here we show the lifetime of a hypothetical sys- 
tem with Sun-Earth-Moon parameters with varying plan- 
etary mass from 2 to 25 Mq. Each dot represents the 
result from numerical solutions. The curve is generated 
by equation (22 1, and equation (24). In this region, the 



lifetime of the moon increases exponentially as the mass 
of the planet increases. The border between Type II and 
Type II is M p = 17.2M®. The border between different 
types depends on the initial conditions and Q p . 



much longer than that of the Sun. The lifetimes of the 
0.3 Mq and the Sun are the order of 100 Gyr and 10 Gyr, 
respectively. Like the 'ice line' with respect to planet for- 
mation, we define the 'moon stability line' as the location 
beyond which a moon is stable for the life of the stellar 
system. Therefore no such primordial moons can presently 
exist inside the moon stability line, though moons are pos- 
sible outside it. Each point represents a moon stability 
line. In each case, we assumed the planet-moon synchro- 
nized state, i.e. n m (0) = £l p (0) as the initial condition, and 
the moon almost reached the outer most stable radius, i.e. 
n m (Tl) w n crit (Tl) and n m (Tl) > n erit (Tl)(Figjl6|). Be- 



cause the moon did not reach the outer most stable radius, 
after that the system continued to evolve back inwards to- 
wards the planet. 

It is worth noting that the moon stability lines de- 
pend on Q. In FigjMj we adopt Q=12 for rocky planets, 
Q=10 4 ~ 10 5 for ice and gas giants. Increasing Q increases 
the tidal evolutionary timescales. 

Overall, the moon stability line moves inward for mas- 
sive planets, for less massive parent stars, and for younger 
systems. In other words, moons are more stable when the 
planet/moon systems are further from the parent star, the 
planets are heavier, or the parent stars are lighter. This re- 
sult can be explained by the size of Hill radius. The planet 
has a larger Hill radius for larger planet mass, smaller stel- 
lar mass, or larger planetary semi-major axis. In general, 
the moon has longer lifetime for the larger Hill radius of 
the planet. 

Gas giant systems For gas giants, the moon stability 
line moves inward as the mass of the planet increases. In 
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Fig. 14. — This graph shows the location required for the planet/moon system to have 1 Gyr, 2 Gyr, up to 10 Gyr lunar 
lifetimes. The blue lines are for 1.0 Mq stars. The red lines are for 0.3 Mq stars. The pictures show which planet-moon 
pair the system has. The size of the pictures do not accurately depict the size of the planet and the moon. The top solid 
lines, both red and blue, are for 10 Gyr lunar lifetimes. The second and third dashed-lines are for 5 Gyr and 1 Gyr in both 
red and blue lines. We used the planet-moon synchronized state as the initial condition for each case, i.e n m (0) = f2 p (0). 
Each star-planet- moon system has ten dots. From the bottom to the top, the dots represent 1 Gyr, 2 Gyr, up to 10 
Gyr lunar lifetimes. For the Earth, we used a k^p of 0.299, a Q of 12, and the moment inertia constant, a of 0.33. For an 
8-Earth-mass-planet, ki v is 0.299, Q is 12, a is 0.33, and the radius of the planet is 1.8 times the Earth radius. For Neptune, 
&2p is 0.4, Q is 10 4 , and a is 0.23. For Saturn, ki p is 0.35, Q is 1.8 x 10 4 , and a is 0.21. For Jupiter and the 10 Jupiter 
mass planet, ki v is 0.5, Q is 10 5 , and a is 0.254. We assume that the 10- Jupiter-mass-planet has the same radius as Jupiter. 
Note that the moon stability lines depend on Q. We adopt Q=12 for rocky planets, Q=10 4 ~ 10 5 for ice and gas giants. 
Increasing Q increases the tidal evolutionary timescales. 



other words, moons are more stable when gas giants are 



heavier. Barnes & O'Brien (2002) studied this type of sys- 
tem. They concluded that smaller moons are more stable 
around gas giant planets. With their high mass ratios, 
from these results we agree that smaller moons are more 
stable around heavier gas giants. 



Ice giant systems In Fig 15 top right, we use &2 P = 0.4, 
a = 0.23, Q p — 10 4 which are the parameters of Neptune, 
and use the Earth as a moon. In this case, the moon 
stability line does not move much even the mass of the 
planet increases. 



Rocky planet systems Compared to ice and gas gi- 
ant planets, rocky planets are small and light. Therefore, 
the mass of the moon is one of the factors that moves the 



moon stability line when compared to the giants. Look 
at the Earth/Moon system and the Earth/Mars system 
in Fig fill The only difference between these systems is 
the mass of the moon. As you can see, the moon stabil- 
ity line moves inward as the mass of the moon increases. 
Look at the Earth/Mars system and an 8- Earth-mass- 
planet/ Venus system in Fig 14 The differences between 
these systems are the masses of the planet and moon. But 
the mass ratio between planet and moon is about 10 to 1, 
which is the maximum ratio of planet and moon for which 
our formulation is valid, in both cases. As you can see, 
the moon stability line moves inward as the mass of the 
planet and moon increase to a ratio of 10 to 1. FigfTB"] top 
left shows the moon stability lines for rocky planets with 
our Moon mass moon. We use &2 P = 0.299, a = 1/3, and 
Q p = 100. The moon stability line moves inward as the 
mass of the planet increases except below the 1 Mgj . 
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Fig. 15. — The top left graph shows the stability line for rocky planets with our Moon mass moon. We use k 2p of 0.299, a 
of 1/3 and Q p of 100. We find the planetary radius from Fortney et al. (2007) with the same composition rate of the Earth. 
The top right graph shows the stability line f or ice giant planets w ith the Earth mass moon. We use k 2p of 0.4, a of 0.23 and 
Q p of 10 4 . We find the planetary radius from Fortney et al. (2007) with the same composition rate of Neptune. The bottom 



left graph shows the stability line for gas giant planets with the Earth mass moon. We use k 2p of 0.5, a of 0.254 and Q p of 
10 5 . We assume that the planetary radius is one Jupiter radius. 



The planet /moon system is preferred to be closer to the 
parent star to detect an extrasolar moon because we can 
make many observations of the planet and moon transiting 
. However, our result shows that the moon stability line 
moves inward for a younger system. If the planet/moon 
system to the parent star is close, we may find the planet 
but it's moon has already gone. If the planet is far from the 
parent star, it's moon may exist but the planet is hard to 
detect. We expect that the semi-major axis of the planet 
around which the first extramoon of a G-type star is 0.4- 
0.6 AU because the lifetime of the moon is more than 10 
Gyr in most cases and we can observe the transiting planet 
two to four times in a year. For M-type star, we expect 
that the planet/moon system locate 0.2-0.4 AU because 
the lifetime of the moon is more than 10 Gyr in most cases 
and we can observe the transiting planet three to six times 
in a year. 

7. CONCLUSION 

We derive analytical expressions for determining the 
lifetime of hypothetical moons in star-planet-moon sys- 
tems. Our solutions allow us to find the type of system 



and the lifetime of the moon without the need to numer- 
ically solve a system of differential equations. We first 
determine whether the moon remains within the planet's 
outermost stable orbit. If not, the moon is lost and the 
system is Type IV. If the moon remains in orbit, there are 
three possible outcomes: Types I, II, and III. In Type I, 
the planet is tidally locked with the moon. In Type II, the 
planet is tidally locked first with the star, and later with 
the moon. In Type III, the planet is not tidally locked with 
the moon. The type of system depends on characteristics 
of the star, planet, and moon (mass, radius, Love number 
Q p , etc.) as well as the initial conditions of the planet and 
the moon. 

Once we determine the system type, we can calculate 
the lifetime of the moon. To find the type of system and 
the lifetime of the moon, we need Tl, which is the time 
when the spin angular velocity of the planet is equal to the 
angular velocity of the moon; See Fig[l] Fig(3] and Fig|4] 
We should use a numerical method to find Tl. 



Our results are extension of Ward & Reid (1973) and 



Barnes & O'Brien (2002). At the range that they con- 



sidered, our results agree to their results. |Ward fc Reid| 
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Fig. 12. — This graph shows the lifetime of a hypotheti- 
cal system with Sun-Earth-Moon parameters with varying 
moon mass from 0.01 to 2 Mj. When the mass of the 
moon is small, the system is Type III. As the mass of the 
moon increases, the system becomes Type II, Type I Case 
2, and then Type I Case 1. There is a minimum lifetime 
of 4.24 x 10 10 years at 0.307 Mp. This minimum arises 
because the mass of the moon has different effects in Type 
I and III. In Type III, the heavier moon has a shorter 
lifetime due to faster tidal evolution. While in Type I, the 
heavier moon has the longer lifetime because the moon has 
grater orbital angular momentum. When these effects can- 
cel each other out, there is a minimum value. The borders 
between the Type III and Type II, Type II and Case 2, and 
Case 2 and Case 1 are M m = 0.181M 5 , M m = 0.842M 
and M m = 0.972Mp, respectively. The border between 
different types depends on the initial conditions and Q p . 



( 1973 ) considered Type III without critical mean motion. 



In this case, the planet will lose its moon only if the moon 



collides with the planet. Barnes & O'Brien (2002) consid- 



ered Type III with critical mean motion. In this case, the 
moon may either hit the planet or escape from it. In both 
cases, the planet and moon are asynchronous. 



Barnes & O'Brien ( 2002 ) concluded that the heavier the 
moon, the shorter the lifetime of the moon. Because they 
considered only systems of Type III, this result agree to 



our result (Fig 12). On the other hand, the heavier the 



moon, the longer the lifetime of the moon for Type I and 
II. 

Our Moon stabilizes Earth obliquity - a key reason for 
the development of life on Earth ( |Ward fc Brownlee|2000 1 . 
Stable obliquity in its star's habitable zone may be neces- 
sary for a planet to support life. An extrasolar moon of 
sufficient mass could stabilize the obliquity of an Earth-size 
extrasolar planet. However even if a planet has a relatively 
large moon like the Earth does, the planetary obliquity 
may not be stable in some cases such as the moon is far 
from the planet, the planet is close to the star, there is 



a Jupiter-size-planet close enough to the planet, etc (Lis- 
sauer et al.||2012 ) . On the other hand, Mars has relatively 
small satellites, and its obliquity changes chaotically, fluc- 



tuating on a 100,000-year timescale (Laskar & Robutel 



10.5 



2 4 

Mass of the moon (M Moo]1 ) 



Fig. 13. — This graph shows that the lifetime of a hypo- 
thetical system using Sun-Earth-Moon present parameter 



with (this work) and without (Barnes & O'Brien (2002)) 



lunar tides in log scale. The black thin line represents the 
lifetime of the moon including lunar tides. The light blue 
thick line is the lifetime of the moon not including lunar 
tides. In the Type III region, both results agree very well. 
However, in the Type I region, the necessity of incorporat- 
ing lunar tides' effect on the planet's rotation, as we have 
introduced in this paper, becomes clear. 



1993). Having a relatively large moon is not enough in 



and of itself to provide a sufficient condition for an extra- 
solar planet to stabilize its obliquity meaning support life. 
Hence, our results give a condition needed to support life 
on a planet in the habitable zone. 

Suppose we find a Jupiter-size planet in the habitable 
zone. This planet may have an Earth-sized moon. If the 
lifetime of that extrasolar moon is equal to or greater than 
the age of Earth, then the moon may support life. Hence, 
our results gives a condition needed for potentially habit- 
able moons. 

In the third example, we show the moon stability lines 
for 1 to 10 Gyr applied to types of planet/moon systems. 
We define the 'moon stability line' to be the location be- 
yond which a moon is stable for the life of the stellar sys- 
tem. In general, the moon stability line moves inward for 
more massive planet, for a less massive parent star, and for 
younger systems. In other words, moons are more stable 
when the planet/moon systems are further from the par- 
ent star, the planets are heavier, or the parent stars are 
lighter. We expect that the semi-major axis of the planet 
for the first extramoon of a G-type star will be 0.4-0.6 AU 
and for an M-type star 0.2-0.4 AU. 

This lays the ground work for the tidal evolution of a 
star-planet-moon system and makes it possible to classify 
star-planet-moon systems and providing useful estimates 
of the lifetime of a moon. 

In some cases, we may not necessarily be able to ac- 
curately predict the long-term survival of the moon. The 
value of Q p , the specific dissipation function of the planet, 
is assumed to be constant in time. However, Q p is not 
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Fig. 16. — This is the first part of a typical graph for 
the third example. The initial condition is the planet- 
moon synchronized state, i.e n m (0) = f2 p (0). The moon 
almost reached the outer most stable radius, i.e.n m (Tl) w 
n crit (Tl) and n m (Tl) > n crit (Tl). The system continues 
to evolve back inwards towards the planet. 



known theoretically, and may depend on the planetary in- 
ternal structure. For the sake of simplicity, we considered 
a star-planet-moon system with a single planet and a sin- 
gle moon. But we do not consider any interactions be- 
tween the star and the moon. This deficiency may be ad- 
dressed in future work. Gravitational perturbations caused 
by other planets or moons may be significant. For close- 
in planets, the stellar gravitational perturbations of the 
moon's orbit are important ( Cassidy et al. 2009 ) . In these 
situations, our method may not predict the lifetime of the 
moon accurately. Despite its shortcomings, our approach 
provides an important step toward understanding the tidal 
evolution and longevity of extrasolar moons, and will form 
both a basis for future theoretical investigations and direc- 
tion for future searches to detect extrasolar moons. 



Appendix 
Type I Solution 

Suppose that the initial conditions n m (0), f2 p (0), and 
rip(0) are known. From these initial conditions, we can 
calculate n rn (t), Q p (t), and n p (t) by solving the equations 
(10) numerically (Fig[I]). 



For < t < Tl, we can use equation ( 13 1 with sgn(£l p — 
n m ) = sgn(f2 p — n p ) — 1 because the planet is not tidally 

> 
1 



locked with either the star nor the moon, and £l p > n. 



For Tl < t < T, we can use (14 1 with sgn(f2 p — n p ) 



because the planet is tidally locked with the moon, and 



^* Tip 



Define a function n m (t) such that n m (t) satisfies the 
following equation for < t < T; 
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The planet will lose the moon. 
Use equation (26). 



YES 



The system is Type I Case 1. 
Use equation (21). 



YES 



The system is Type I Case 2. 
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Fig. 17. — Flow-chart for calculating moon lifetimes in a 
star-planet-moon system. First, check the type of system. 
Then, calculate the lifetime of the moon. 



where Tq 



M m (GM p ) 2/3 

«m (0) 



aR 2 p M p n p (0) 



M p (GM 3 f/ 3 . c 
«p /3 (0) 



the initial angular momentum of the system. In other 
words, n m (t) represents situation in which the planet and 
the moon are tidally locked from beginning to end. Be- 
cause n m (t) and n m {t) are the same for Tl < t < T, we 
can calculate the maximum lifetime of the moon if we know 



the domain of n m (t) (Fig 18) 
For given t, set x — 



and define 



/(*) 



_ (GM m )(GM p ) 2 / 3 



,1/3 



+aR 2 p (GM p )x- 



(GM P )(GM S ) 2 / 3 
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The condition that f(x) has at least one zero, i.e. equation 
( 14b ) has a real solution, is 



(38) 



Since we know n p (t), we can plug in equation (14a) and 
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But in this equation, we measure the time t from Tl. Mea- 
suring the time t from zero, for Tl < t < T2, 



n rn (t) 
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Fig. 18. — The thick dashed line is n m (t). The thick 
light blue line is n m {t). As you can see, n m (t) and 
n m have the same maximum time. We use the present 
data of our Sun-Earth-Moon system. The initial condi- 
tions are n m (0)=84 rad/year, fi p (0) — 730ir rad/year, 
and n p (0) = 2tt rad/year. The new initial conditions 
are n m (0)=17 p (0) = 48.5524 rad/year and n p (0) = 2ir 
rad/year. 
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n m > from Fig. [3j Hence, 
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Combine equation (41), equation (42), and equation 



@. Solve for T2. 

T2=T1 + , 2 Q » < GM P' 8/3 
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(44) 



solve for i: 



Next, we will find n pc and T3. At t = Tl + T2, we can 
see that f2(t) = n p (t) = n pc from Fig|3] By the conserva- 
tion of angular momentum and equation (|41|) , 



t < 
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Hence, the lifetime of the moons for Type I, T, is 
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(45) 

The second term, aR 2 {GMp)n pcl is the spin angular mo- 
mentum of the planet. We will ignore this term to approx- 
imate n pc because we know the spin angular momentum is 
small compared to the total angular momentum. 

( 40 ) Hence, 



Type II Solution 

For Type II, there are three stages (Fig{3j). We start by 
finding T2. Assume n m (Tl) and n p (Tl) are known. At 
the end of the planet-star synchronized state, the torque 
due to the star is equal to the torque due to the planet. 
Hence, r p _ s (Tl + T2) = r p _ m (Tl + T2). From equation 
([I]), equation ^ and Kepler's Law, 



f ( GM m ) 7/ 6 ( GM P ) 1 /2 + (GM P )(GM S ) 2 / 3 
I GL 



(46) 



When we use O p (0) = n p (0) = n pc and n m (0) = 
n m (Tl + T2) as a initial condition, we can calculate T3 
by using the formula for Type I. Hence, 
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where n pc = n p (Tl + T2) (Fig. [3}. 
For Stage 2, we can see that f2 p 
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From equation (44), equation (46), and the equation 



(47), the total life time, T, is 



Hence, by the equation (15a) with n m (0) = n m (Tl) 
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Fig. 19. — This is the graph of n m (t) in Type III. In this 



graph, A 



_ 39 k 2p R 'p GM m 
2 Q p (GM p )8/3- 



of the torque due to the star at t = Tl (Figjl]), 

\r p - m (Tl)\ > \r p - s (Tl)\. (52) 

By equation 0, equation ([2|, and Kepler's Law, equation 
( 52 1 implies that 



(53) 



For the time period < t < Tl, the planet is not tidally 
locked with either the star or the moon. We can use equa- 
tions (13 1 with sgn(f2 p — n m )=l and sgn(r2 p — n p )=l. 



Combine equation (53), equation (13a), and equation 



(13b), then solve for Tl. We obtain 
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Type III Solution 

In Type III, the graph of n m (t) is comprised of two 
parts. The first part, n^t), is from < t < Tl and the 
second part, n~ (t) is from Tl < t < T. sgn(£l p — n rn ) is 1 
for < t < Tl and -1 for Tl < t < T because the planet 
and the moon are not in a synchronized state. To have 
n+(i), we can use equation (13) directly. To have n~(t), 



we set n m (0) = n m (Tl) and t = t — Tl because we need to 
shift the graph Tl in a positive direction. From equation 
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for Tl< t < T. 
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We can get the same result if we set the inside of the paren- 



thesis of equation (50) equal to 0, and then solve for t. 



Condition for Type I Case 1 

The condition for Type I Case 1 is that the magnitude of 
the torque due to the moon is greater than the magnitude 



Condition for Type I Case 2 

Assume n p {Tl) and n m (Tl) are known. Let be the 
time from Tl, when the magnitudes of two torques are 
equal (Fig[7|. 

The condition for Type I Case 2 is 

fip(t.) > n p {U) (55) 

where i* satisfies 

\T p - m {t*)\ = \r p - s (U)\. (56) 



Equation ( 56 ) implies that 
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We did a similar calculation to arrive at equation ( 53 1 
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By conservation of angular momentum and equation 
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After Tl, the planet is not tidally locked with either the 



star or the moon for a time. We can use equations ( 13 ) 



with sgn(51 p — n m ) = — 1 and sgn(il p — n p ) = 1. We use 
equation (13a) and equation (13b) with initial conditions 



n>m(0) = ^m(Tl) and n p (0) = n p (Tl) , which are 



n m (t) 
and 

n p (t) = 



39 k 2p R 5 p GM m 
Y Q p (GM p )8/3 



-3/13 



t + n-^{Tl) 



39 k 2p R 



2 Q p (GM p )(GM s )Vz p 



t + n-^iTl) 



(59) 

-3/13 

(60) 



17 



Plug in equation ( 59 1 and equation ( 60 1 into the equa- 



tion (|57|), and solve for t*. Then, we have 
t 
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By using this £*, we have 
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Applying equation (55), we have 

ai 6 12 X 4 - GLo* 3 + -6 3 < 0. 



Conditions for Type II and III 
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From Fig|8]and Fig[9j the condition for Type II is T5 > 
T6 and the condition for Type III is T5 < T6. 

For Stage 3, we know 
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At t = T6, n m (t) = oo. Hence, 
2 g p (GM// 3 _ 
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To find T5, we use the same formula for Type I, which 
is equation Q. We set n p (0) = n p (Tl + T2). T5 is 
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From equation (|46|), we also know that 
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The condition for Type II, T5 > T6, implies 
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Similarly, the condition for Type III, T5 < T6, implies 
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1/4 



13 



(71) 



((GM m ) 7 / 6 + (GM p ) 1 /2( G M s ) 2 / 3 
< 3 39 / 4 (GM p ) 6 (GM s ) 8 (GL ) 13 . 
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